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The role of Klein factors is investigated for the bosonized Hamiltonian of the dimerized Hubbard model. 
Contrary to previous approaches we take into account their number changing property, i.e. we do not replace 
them by Majorana fermions. We show how to treat Klein factors in the framework of the self-consistent 
harmonic approximation, both for finite systems and in the thermodynamic limit. 



1 Introduction 

The foundations of bosonization were laid more than 50 years ago in a seminal paper by Tomonaga |fT|. 
During the following decades the method was worked out and successfully applied to one-dimensional 
electron and spin systems |2||3]E1- Despite its long history there are still some subtle points in the bosoniza- 
tion formalism which are not taken into consideration in the majority of the literature. One of these is the 
proper treatment of the so-called Klein factors which have to be introduced in order to preserve the an- 
ticommuting property of the fermionic fields during the bosonization procedure. The role of the Klein 
factors deserves particular attention when nonlinear perturbations arising e.g. from impurity scattering or 
lattice modulations are to be considered in finite systems. 

In this paper we demonstrate how to handle the Klein factors in a systematic way, both in the thermo- 
dynamic limit and for finite systems. As a prototypical model we study the one-dimensional dimerized 
Hubbard model where the hopping is periodically modulated due to the Peierls distortion of the lattice. We 
extend the self-consistent harmonic approximation |5||6|E1 by treating bosonic fields and Klein factors on 
equal footing [8J. As an application we use the formalism to calculate spin and charge gaps of this model. 

2 Bosonization and Klein factors 

We consider the dimerized Hubbard model 

which differs from the ordinary one-dimensional Hubbard model by a periodic modulation of the hopping 
described by the dimerization parameter u. This modulation is relevant in one dimension due to the cou- 
pling between the lattice and the electronic degrees of freedom. A finite u corresponds to a periodic lattice 
distortion known as Peierls instability. We study the case of half-filling, and have chosen the modulation 
accordingly to be of the form {—\yu. The quantity u is considered to be a parameter of the model, 
creates an electron with spin direction a =t, i at site i, and t and U are the hopping matrix element and 
the on-site Hubbard interaction, respectively. 

In the following we list the main steps which have to be performed in order to bosonize the Hamiltonian 
Q. For more details we refer the reader to the reviews j3||4j|6|. First we represent H in momentum space: 
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H = ^efcC+^Cfc^ + to^(e*'=c+,Cfe+^_^ + h.c.) + j^Y^ c++^,tCfe|C+_^_^Cfe,^ (2) 

/c,(T fe,cr k,k' ,q 

where e/j = —2tcosk, and N is the number of lattice sites. We hnearize the spectrum around the two 
Fermi points ikp, kp = 7r/2, and introduce left and right moving fermions labeled by L and R, re- 
spectively. This allows us to sort the various scattering processes according to their initial and final states 
("g-ology"). Then we define bosonic operators 6+^ which are related to fermionic particle-hole excitations 
via ["^L/Rcrik) = Cq:(fep+fc),cr ctc; Cj.^, c^^ are the standard Fermion operators] 

bta = -^J2mk + q)-^aik) (9>0) (3) 



where q = {2Tr/N)nq and a = i?t,Lt,iiJ,,LJ,. Returning to real space we define the bosonic fields 



9>0 V*^? 9>0 V"9 

which are related to the fermionic field operators via the bosonization identities 

^La{x) = ^Fi<,e-''^i-(^)e-*^^<'(^)e-2'''^^<'^/-^ (5) 
yL 

= ^ p^^(,-i(vt^M+(pLaM)^-2wiNLaX/L 

■v/27ra 
\/2na 

where Na counts the particle number with respect to the filled Fermi sea, and L is the length of the system. ' 
The Klein factors F+ (_F„) are unitary operators that commute with the bosonic fields. They change the 
number of the fermion species a by ±1, a change which cannot be achieved by any combination of the 
bosonic field operators. In order to ensure the correct anticommutation relations for ^a{x), the Klein 
factors have to fulfill the following relations: 

{F^,Fp} = {F+,F+} = 0, {F+,Fp} = 26c,p, [N^,Fp] = -d^^pF^ . (8) 

Here [., .] denotes the commutator, and {., .} the anticommutator. In a last step one combines ipa and 
to new fields (j)c^s and ^c,s. 

{(fL^ ± ifiLi + vm ± V>Rl + ii-c-) (9) 



2V2 

dc,s = -^-^{(PL^ ±(PLI - fRT T fRl +ii-C.) (10) 

and introduces TV,; .5 = ± Ni and Jc,s = ± Jj^ where N^r = N^a + ^Ra and Ja = -^Lu — Nr^. 
Combining everything we obtain 

H = Ho + Hi+H2 + H3 (11) 

where Hq is the Luttinger Hamiltonian: 

^0= E ['^■■{-id.M' + Va9.idJc.A:+{^ E {-K + VagaA I (12) 



^ We have chosen the lattice constant to be unity, i.e. L = N. Nevertheless we retain L and N for easy reference. 
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Hi the Umklapp contribution: 



UL 



^'^^Jo '^'^ ^^«T^^i^ii^iTe-^'^'^= + h.c.} , U = ; (13) 

H2 the backscattering contribution: 

H2 = U dx {F+F+^FR^FL^e-^^^^^ + h.c.} ; (14) 
Jo 

and Hs the dimerization contribution {ii = tu/TTo): 

H3 = u dx {tF+FLiie-'^^^^+^^^ + iF+^FLie"^^*^-'^^^ + h.c.} . (15) 
Jo 

In Eqs. ( I13> - M5\ we have used A^c = = at half filling. The parameter a is a short-distance cutoff 
of the order of the lattice spacing, i.e. of order one. The Luttinger parameters gc.s and the charge and spin 
velocities Vc.s can either be calculated perturbatively or from the Bethe ansatz solution of the Hubbard 
model I9I IIOI . In the following we focus on the role of the Klein factors. In the Uterature it is common 
practice either to ignore them or to replace them by Majorana fermions f4l II II . It is argued that the 
latter approach - which neglects the number changing property of the Klein factors - should be justified 
in the thermodynamic limit. In the following we aim to present a more rigorous approach. Due to the 
conservation of charge and spin all combinations of Klein factors appearing in the Hamiltonian (II 1> can 
be expressed in terms of the operators — F^^Fl] and Ai — F^^Fli plus their hermitean conjugates. 
In particular, the four-fermion terms arising from Umklapp and backscattering read 

Pr-\-PriPliFl^ = F^-^FltF^^Fli = A^Ai , (16) 

Since the Klein factors are unitary, F^Fa — F^F^ — 1, it is easy to show that 

[A^,A+] = [A^,A+]^Q , (18) 

and we may choose a basis where A„ and A'^ are both diagonal. From yl^ A| — A'^ A^ — 1 one concludes 
that the eigenvalues of Aa are pure phase factors, i.e. 

A^\k^,kC) = e^'^Mfct-fci) A+\k^,kC) - e-"=T|fcT,fc;) (19) 

A^\k^,k{) - e^^Mfct-fci) = e-"=MfcT,fci) (20) 

with < fccr < 2n. The terms ^ s appearing in Hq do not commute with the Klein factors; however 
it appears reasonable to neglect them in the thermodynamic limit L — > 00. We will come back to this 
question in Sec.|3 We thus replace the Klein factors in Hi , H2 and H^ by their eigenvalues, and obtain 

Hi = U dx {e^('=T+fei)e-*2V20e ^ y ^21) 
Jo 



i-L 

H2 = U dx {e^C^T-fei) 6-^2^2^= + h.c.} 

Hs^ii dx {ieik,^-iV2{^.+4,.) ^ jeifcig-*V2(</.c-0.) _^ ^i.c.} . 
Jo 



(22) 



(23) 
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As a result the Hamiltonian of the dimerized Hubbard model separates into different sectors of purely 
bosonic Hamiltonians which are labeled by fc| and . Note that when replacing Klein factors by Majorana 
fermions [4 | one obtains only the eigenvalues ±i for the two-fermion terms and ±1 for the four-fermion 
terms, i.e. continuity is lost. Shifting the field operators according to 0c, s 4>c,s + [k] ± ki)/2^/2, the 
phase factors can be absorbed, with the result 

H ^ Ho + 2U dx (cos 2 V20c + cos 2 V2(/)s) + I dx sin V2(/)c cos V20s . (24) 
Jo Jo 

In this sine-Gordon-like Hamiltonian the operator constraint [0c s]q=o = 0, see Eq. (0}, has to be replaced 
by 

dx <i)c,s^ ^^^{k]±ki) . (25) 

3 The self-consistent harmonic approximation 

In order to study the charge and spin gaps in the dimerized Hubbard model we use the self-consistent 
harmonic approximation (SCHA) in which the exponentials of field operators appearing in Mil - ( I23> are 
replaced by quadratic forms. We introduce the trial Hamiltonian 

Htr = + v^g^{d,9^)^ + ^0^1 (26) 

which provides us with a variational estimate for the ground state energy 

E = {Ha)tr + {Hi)tr + {H2)tr+{H3)t, , (27) 
E Etr sr-^ Al{(pl)tr . , 

T = ^^^+^(^=^,^1). (28) 

Here the expectation value is with respect to the ground state of Htr, Eti is its ground state energy, and 

E{k^,ki) = 2Bicos(fc| + fcx) + 2B2Cos(fc| - fcx) - 2B(sinfct +sinfcx) (29) 

with 

B^ = i} e-4<'^c>tr ^ B2 = U e-^^-^'^- , B = u e-<*c>tre-(*^>t. . (30) 
Minimizing E with respect to Ac and A^ yields the gap equations 

^ . -4B,^-b'J^ (31) 

27rWc5c oBi oB 

^ . -4a||^-B^ (32) 

zTTVsgs 0B2 oB 

where Eq is the minimum of E{k^, k^) (see Table 1). In order to solve these equations analytically we 
consider the case U > where gs — I and gc < 1, and restrict ourselves to the limit of small dimerization 
u. Since iJtr is quadratic in the bosonic fields it is straightforward to calculate 



a2\ 



9c , Ao 



/tr 



In-^ (33) 



2 A 



^Ihr = ^In^ (34) 
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range 






Eo{Bi, B2, B) 


< B < 2^2 


arcsin 


TT — arcsin 


-2Bi - 2B2 - %^ 


2^2 < B 


2 


TT 

2 


-AB - 2Bi + 2B2 



Table 1 Minimum £0 of i5(fc| , fc| ) - see Eq. <29t - whicii is used in the gap equations <31> and <32t . The first line 
is applicable for u = Q (Hubbard model without dimerization) whereas the second line applies for it > 0. 



where Aq is a cutoff-dependent energy scale of the order of the bandwidth, and Ac,s < Ao is assumed. 
For nonzero dimerization u > the solutions of the gap equations lie in the range B > 2B2; thus the 
second line of Table 1 with dEo/dBi = —2, dEo/dB2 = 2 and dEo/dB = —4 has to be used in Eqs. 
J3U . ( I32t . For u ^ the spin gap vanishes while the charge gap approaches a constant according to 1121 

Ae(M)-Ae(0) CX U^/^ (35) 

As{u) CX u^/^ (36) 

with cutoff-dependent prefactors. The exponent 2/3 that characterizes the spin gap is in accordance with 
the corresponding exponent of the dimerized antiferromagnetic Heisenberg chain up to a logarithmic cor- 
rection in the prefactor 1131 1141 . Since the Heisenberg model corresponds to the U ^ oo limit of the 
half-filled Hubbard model, this indicates that as far as the exponent is concerned the SCHA result ( I36l l is 
exact and persists even in the strong-coupling regime U/t 3> 1. For u > Uco the behavior of the gaps is 
changed to 1151 

Ac(u) w A,iu) CX u2/(3-9c) (37) 

where the crossover value Uco is defined by As{uco) — Ac(0). In Fig. 1 we show Ac{u) — Ac(0) and 
As{u) as a function of u for U/t = 2 as obtained from the numerical solution of the gap equations. For 
comparison, the analytical results ( I35t and ( I36t are also given. 



4 Finite systems 

For a finite system of length L it is not possible to simply replace the Klein factors by their eigenvalues 
since the terms ~ ^ in the Luttinger Hamiltonian il2\ do not commute with the F's. However one may 
decouple the Klein factors from the bosonic fields using a variational ansatz. To this end we introduce the 
"Klein Hamiltonian" 

Hu- = iB{F^^FL] + F^iFli) + BiF^^F^^FliFl] + B2F^^F^^FriFl] + h.c. 

+^{vcgcJl + VsgsJl) (38) 

where Bi , B2 and B are variational parameters to be determined self-consistently. H^^ is of the form of a 
tight-binding Hamiltonian for a particle moving on a 2d lattice in a harmonic potential. A class of similar 
Hamiltonians for general potentials in Id has been studied in 1161 . The choice of the trial Hamiltonian is 
equivalent to a decoupling of the non-linearities according to 

{F+^F+^FR^FL^),,e'^^^^^' + F+ F+ i^Ki^^LT (e-'2^*=)tr , (39) 
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Fig. 1 Charge gap Ac and spin gap As (in units of t) of the dimerized Hubbard model for U/t = 2 obtained within 
the SCHA. [We used tic, s ~ vf\/1± U/ttvf, gs = 1, and gc = + U/tvvp; note that i^f = 2f]. The straight 

lines are the analytic resuhs <35> and <36> valid for u < Uco- Note that Ac(0) is subtracted from the charge gap in 
order to highlight the power law behavior. 



and analogously for the Umklapp and dimerization term. This means that instead of replacing the products 
of Klein factors by their eigenvalues as in the thermodynamic limit, we now have to replace them by their 
expectation values with respect to the ground state of H^. In both cases one ends up with a sine-Gordon 
type model like the one explicitly given in Eq. \2A\ . It is conceivable that for a finite-size system the 
eigenvalues are not identical to the expectation values, i.e. different sine-Gordon models arise. 

The question is then whether identical sine-Gordon models arise at least in the thermodynamic limit. 
We can answer this question within the framework of the SCHA. Here the introduction of the "Klein 
Hamiltonian" amounts to replacing the quantity Eo{Bi, B2, B) which enters the gap equations (and which 
is explicitly given in Table 1) by the ground state energy of H^. Apart from this modification Eqs. J30> 
- ( I32> remain unchanged. Consider first a system with both a spin and a charge gap. For large systems 
the parameters B, Bi, and B2 become size-independent, with the consequence that the kinetic energy in 
dominates the confining potential. The ground state energy of is then given by the minimum of 
E{k^ , k^), see Eq. j29> . and expectation values of the Klein factors are equal to their eigenvalues. 

In the absence of dimerization, i.e. for the "pure" Hubbard model, we encounter a different situation 
since the spin gap vanishes. For large system size one finds B = 0, Bi = const, and B2 oc U{a/L)^. 
Thus the confining potential in the spin sector dominates the corresponding kinetic energy (proportional 
to B2) in the thermodynamic limit. Hence the ground state of is the eigenstate of the current operator 
with Jg = 0, and the expectation value of the operator F^^F^^Fh^Fj^t^ in this state is zero, in contrast to 
the eigenvalues which are exp(i/c| — ik^), see Eqs. (I16> - ilOi . 



5 Conclusions 

We studied the role of Klein factors for the bosonized Hamiltonian of the dimerized Hubbard model. Since 
the Klein operators do not commute with the total spin and charge currents Jc s they cannot acquire a 
definite value at the same time as the current operators. 

The ground state of the gapped system is a superposition of many spin and charge states. In this situation 
it is justified to choose a fixed phase for the Klein operators. The bosonized Hamiltonian is then the 
conventional sine-Gordon-like Hamiltonian, cf. Eq. (I24> . In an ungapped system where the non-linearities 
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introduced by backscattering are (marginally) irrelevant operators, the ground state has a well defined 
current. In this case the phase of the Klein operators is undetermined and Eq. (I24t has no justification. 

We have extended the self-consistent harmonic approximation in such a way that Klein factors can be 
handled systematically, and we worked out the theory for the Hubbard model. In a previous paper (8| we 
applied the SCHA to spinless fermions with nearest-neighbor interaction and dimerization, and studied in 
detail finite-size effects. In the spinless case the Klein Hamiltonian can be mapped onto a Mathieu 
equation which can be solved analytically in certain limits. It turns out that finite-size corrections to the 
gap equation are not important as long as the dimerization gap is larger than the finite-size gap. In addition, 
the finite-size formalism allows to calculate the size-dependence of the Drude weight in the gapped regime 

m 
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